Muon Fluxes from Dark Matter Annihilation 
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We calculate the muon flux from annihilation of the dark matter in the core of the Sun, in 
the core of the Earth and from cosmic diffuse neutrinos produced in dark matter annihilation 
in the halos. We consider model-independent direct neutrino production and secondary neutrino 
production from the decay of taus produced in the annihilation of dark matter. We illustrate how 
muon energy distribution from dark matter annihilation has a very different shape than muon flux 
from atmospheric neutrinos. We consider both the upward muon flux, when muons are created in 
the rock below the detector, and the contained flux when muons are created in the (ice) detector. 
We contrast our results to the ones previously obtained in the literature, illustrating the importance 
of properly treating muon propagation and energy loss. We comment on neutrino flavor dependence 
and their detection. 



I. INTRODUCTION 

The dark matter problem, where more matter is re- 
quired to account for gravitational forces observed on as- 
tronomical objects than it is visible, has persisted for 
more than seven decades [l|. Observations of galactic 
rotation curves Q, orbital velocities of galaxies within 
clusters Q, anisotropics of the cosmic microwave back- 
ground (CMB) [H, distance measurements from Type la 
supernovae (SN) and baryon acoustic oscillations (BAO) 
H, and large scale structure @ all imply the existence 
of cold (non-relativistic) dark matter (CDM) with an 
abundance of 23% of the total density of the Universe 
(CIcdm — 0.233 ± 0.013). In addition, the combination 
of the CMB, SN and BAO data predicts that only 4% of 
the total density of the Universe can be attributed to the 
baryonic matter (flbaryons = 0.0462± 0.0015). Thus, the 
particle content of the CDM can not be explained in the 
context of the standard model. 

In all extensions of the standard model, there are many 
candidates to account for CDM: weakly interacting mas- 
sive particles (WIMPs), axions, superheavy dark matter 
(WIMPZILLAs), and solitons (Q-balls) 0, 0, i, ©] to 
name a few. Among these possibilities, a WIMP of mass 
of order 100 GeV provides a natural explanation for the 
observed density of dark matter today [T(| EH- These 
WIMPs were abundant and in thermal equilibrium in 
the early Universe and then eventually "froze out" due 
to the Hubble expansion. 

An interesting coincidence, independent of the dark 
matter issue, is that the 100 GeV scale is the charac- 
teristic scale of new physics beyond the standard model 
according to naturalness arguments (l2^ . Collider ex- 
periments such as the Large Hadron Collider (LHC) at 
CERN will explore this new scale physics in the near 
future [Til]. The detection and characterization of dark 
matter particles is possible in these LHC searches. How- 
ever, the LHC experiments will not be able to determine 
detailed properties of these particles such as whether they 
are stable and what their couplings are to other particles. 

Apart from the colliders, there are two independent 



but complementary approaches to search for dark matter: 
direct and indirect detection [TJ], including dark matter 
accumulation in the Earth and Sun [Til, O E3] and in 
the galactic center |18| . and subsequent annihilation to 
neutrinos [Tfl, [2(| l2ll l22fl. There are a number of direct 
detection experiments |23l. l24l . [2o1 . [2r| . Direct searches 
can provide valuable data on the dark matter's couplings 
to the standard model. They all look for energy deposi- 
tion via nuclear recoils from WIMP scattering by using 
different target nuclei and detection strategies, and ex- 
pect to observe the same WIMP mass and cross sections. 
Currently, the strongest upper bounds (~ 10 pb) on 
the spin independent WIMP-nucleon cross section of a 
WIMP with mass - 100 GeV come from the XENON 
Dark Matter Search (XENON) (3 and the Cryogenic 
Dark Matter Search (CDMS) (25( experiments. So far, 
contrary to the null results of all other direct searches, 
the DArk MAtter (DAMA) collaboration has observed 
an annual modulation in their data [261 ] which is claimed 
to be due to dark matter particles in the galactic halo. 
Lately, some models have also been proposed to account 
for that modulation signal [27], l28( . 

On the other hand, indirect searches for WIMPs 
through their annihilation (or sometimes decay) into 
standard model degrees of freedom such as positrons, 
antiprotons or 7 ra ys, which has been explored in 
several experiments |29l . [30L l3~il HU, HH and neutri- 
nos with experiments such as Antarctic Muon And 
Neutrino Detector Array (AMANDA) 0, IceCUBE 
[H[, Cubic Kilometer Size (KM3) Neutrino Telescope 
(KM3NeT) (36|. Observations in the recent years 
such as the excess in the positron fraction reported 
by High Energy Antimatter Telescope (HEAT) [30lj . 
the Payload of Antimatter Matter Exploration and 
Light-nuclei Astrophysics (PAMELA) [311, Advanced 
Thin Ionization Calorimeter (ATIC) [32| and Polar 
Patrol Balloon and Balloon borne Electron Telescope 
with Scintillating fibers (PPB-BETS) [H[, an excess 



in microwave emission around the galactic center 37] 
(also called the "WMAP Haze"), a bright 511 keV 
gamma-ray line from the Galactic Center region from 
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INTErnational Gamma Ray Astrophysics Laboratory 
(INTEGRAL) [Hj|, and an excess in the flux of 1-10 
GeV diffuse galactic 7 rays from Energetic Gamma Ray 
Experiment Telescope (EGRET) data [39[ have made 
researchers more excited in their quest for dark matter. 
Theoretical studies of the indirect dark matter detection 
via neutrino signals has recently received a lot of atten- 
tion mull. 

In this paper, our focus is on the muon energy distri- 



butions from v., 



in neutrino telescopes due to an- 



nihilation of WIMPs which are captured in the core of 
the Earth (or the Sun) via gravitational interaction, or 
from annihilation of relic neutrinos [lij. As a result of 
these annihilations, neutrinos are produced at energies of 
the order of the mass of the WIMP and they interact on 
their way to the detector producing an observable muon 
flux. There is an extensive literature on WIMP annihila- 
tion in the Earth's or Sun's core @, EE HH HH] • Here, we 
present a systematic way of calculating this muon flux 
for a few choices of annihilation channel. Our results 
can be used to determine the muon flux as a function of 
energy for a specific dark matter model by summing all 
the contributions from each annihilation mode weighted 
with corresponding branching fractions. We also com- 
pare our muon energy distributions with those obtained 
using other theoretical frameworks widely used in the lit- 
erature [rJ m M, H EH • 

The muon neutrino flux from weakly interacting dark 
matter annihilation (xx annihilation) depends on the 
dark matter capture rate and the dark matter annihi- 
lation rate. In the next section, we review the standard 
evaluation of the (muon) neutrino flux. This is followed 
by a discussion of the theoretical framework of muon sur- 
vival probabilities and the resulting muon flux. Results 
are shown in Section IV, followed by our conclusions in 
Section V. The Appendix includes details of the muon 
neutrino energy distribution from various decay modes 
of fcrmions F in XX ~ * FF where F = v, t, c and b. 



II. NEUTRINO FLUX FROM DARK MATTER 
ANNIHILATION 

The dark matter particles can be captured in the core 
of the Sun or the Earth by interacting with the nuclei 
in the medium. This results in a WIMP density in the 
core that is considerably higher than in the galactic halo. 
The capture rate (C) depends on the composition of the 
medium, the WIMP-nucleus interaction cross sections 
(<7q), the WIMP mass (rn x ), the local dark matter den- 
sity (p x ) and velocity (v) distribution of the WIMPs in 
the halo. After being accumulated in the core of these 
dense objects, the WIMPs annihilate with rate Ta into 
standard model particles which may further decay into 
neutrinos. These neutrinos can reach Earth-based detec- 
tors and create fluxes of charged leptons as a consequence 
of neutrino charged-current (CC) interactions. 

The resulting fluxes depend on how the capture and 



annihilation processes have occurred initially, however, in 
equilibrium these two processes are related: for every two 
WIMPs captured, one annihilation takes place so Ta — 
C/2. This equilibrium condition leads to a maximal flux 
which depends on the capture rate given by 0, EE E3 , 
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The summation in Eq. |T]) is over all species of nuclei in 
the astrophysical object, mjv 4 is the mass of the ith nu- 
clear species with mass fraction fi relative to the Sun (or 
the Earth). The kinematic suppression factor, denoted 
by S(m x /mN i ), for a capture of WIMP of mass m x from 
a nucleus of mass m # 4 is given by 0, EE, E3 



S(x) 



A 



1.5 
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where 



A(x) 



2 (x-1) 2 



< v esc > 



(4) 



(5) 



For the Sun, < v esc >= 1156 km/s and for the Earth, 
< v esc >= 13.2 km/s. We also note that S(x) — > 1 for 
x — ► 1, which means that the capture process is kinemat- 
ically suppressed if m x differs from m N i , and there is no 
suppression if these masses are the same. 

The other quantities in the capture rate expression are 
the form factor suppression Fi{m x ) and the velocity dis- 
tribution function </>j of the ith element. The former one 
is due to the finite size of the nucleus which disrupts the 
coherence in the scattering process, thus, the form factor 
suppression is a negligible effect for capture from scatter- 
ing with hydrogen and helium whereas it becomes impor- 
tant for larger nuclei. The velocity distribution function 
4>i depends on the velocity distribution squared of the 
element averaged over the volume of the astrophysical 
object (< vf >) and is given as @, EE E3]> 



< vf > 



(6) 



If the massive astrophysical object is far from equilib- 
rium, which is most likely the case for the Earth, the 
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annihilation rate is not only dependent on the capture 
rate but also on the annihilation cross section (a) via, 



C 



T A = — tSLnh 2 (t y/CC A ) 



where 



C A = 



< <7V > 



(7) 



(8) 



and Vcs is the effective volume of the core of the Earth or 
Sun, while to is the age of the solar system. It is obvious 
from this relation that the equilibrium condition holds 
only when t y/CCA ~> 1. 

The flux of neutrinos of flavor i from dark matter an- 
nihilation into standard model particles can be written 
as 



dE v 



AttR 2 



5> 



dE v 



(9) 



where (dN v / dE v )p,i is the differential energy spectrum of 
neutrino flavor i from production of particles in channel 
F. In general, this energy spectrum is a function of the 
neutrino energy E v and the energy of the produced par- 
ticle, Ei n . The differential neutrino energy spectra from 
a few dark matter annihilation channels are given in the 
Appendix. The quantity R is the Sun-Earth distance for 
neutrinos produced in the core of the Sun, or the radius 
of the Earth for the neutrinos created in the core of the 
Earth. The sum in Eq. ([9]) is over all annihilation chan- 
nels F weighted with corresponding branching fractions 
B F . 

Neutrinos can be detected via their charged current 
interactions near or in the detector. To avoid the down- 
ward muon background, upward events where neutrinos 
interact with the nucleons in the rock below the detector 
producing muons which then travel up through the detec- 
tor are one category of events. The other are contained 
events, in which the muon neutrinos produce muons in 
the detector ice. In the following sections, we focus on 
evaluating muon energy distribution from interactions of 
the neutrinos produced in the annihilation of the dark 
matter in the core of the Earth, in the core of the Sun 
and cosmic diffuse neutrinos from dark matter annihila- 
tion in the halos. 



III. MUON FLUX 

The muon flux from muon neutrinos from DM annihi- 
lation depends on the flux of muon neutrinos as calcu- 
lated with Eq. ([5]) and attenuation, tau neutrino regen- 
eration and neutrino mixing in transit to the detector. 
For the energies of interest, neutrino mixing for DM an- 
nihilation in the Earth's core is not important. For an- 
nihilations in the Sun, neutrino mixing and tau neutrino 
regeneration may affect the flux of muon neutrinos when 



m x is large. We neglect these effects because we consider 
the value of m x for which only moderate modification of 
the muon neutrino flux is expected [To| . 

For upward events where the muon is produced outside 
the detector, muon energy loss is important. The most 
straightforward evaluation is for contained events, so we 
start with this case. 

In the sections below, we focus on the neutrino in- 
duced flux of muons, produced either in the detector or 
near the detector with muon energy loss included. The 
detector may be modeled by an effective area. For Ice- 
Cube, an effective area for muons, A e ff(E u ,9), can be 
simply parameterized as a function of the muon energy 
at the detector 1431. 



A. Contained Events 

Contained events involve neutrino conversions within 
the detector volume. Denote the muon neutrino flux from 
a source of DM annihilations in the Earth's core or the 
Sun's core at location R from what is effectively a point 
source by 



r{E v ,R) ■ 



dE. 



The probability of the conversion of a neutrino with en- 
ergy E v into a muon with energy E^ over a distance dr 
through CC interactions is given by 



dPcc^drdE, 



Pp- 



dtr Z{E v ,Ep) 
dE u 



(P 



(10) 



Here, p p and p n are the number densities of protons and 
neutrons in the medium, respectively. We assume that 
P P = Pn = ^NaP where Na - 6 x 10 23 is Avogadro's 
number. The differential cross sections da^ n jdE^ are 
the weak scattering cross sections of (anti- ] neutrinos on 
nucleus, which can be approximated by [1 01 ] 



daV 
dE, 



ii 



2m p G F 




(11) 



with a™'P = 0.25,0.15, b^P = 0.06,0.04 and ap p = bP> n , 
b^' p — d p / ' n . The contained event rate, for a detector with 
size £, is 



d(j)p 
dE„ 



R+e 



IR 
+ {V- 



dr 



dE v 



dP, 



cc 



dr dE u dE, 



■{E V ,R) 



(12) 



where the neutrino flux is essentially independent of po- 
sition in the detector given the scale of the Earth. 

The neutrino flux from DM annihilation in the Earth's 
core is not attenuated, to a good approximation, until the 
neutrino interaction length approaches the radius of the 
Earth. This occurs at a neutrino energy of approximately 
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100 TeV. The neutrino flux in Eq. ([12]) is given by Eq. 
©, with R = Re — 6400 km, the radius of the Earth, 



dNu 
dE v 



(13) 



F.u 



The density of the Sun is such that one needs to take 
into account neutrino attenuation due to its charged- 
current interactions as they pass through the Sun. In 
our calculations, we approximate attenuation with the 
exponential decrease in the flux over a distance Sr' in 
the Sun, assuming that the composition of the Sun is 
mostly elemental hydrogen which has mass mu = 0.931 
GeV. 

^-{r' + Sr') = exp(-p(r')a C cSr' /m H )^-{r') (14) 

where the neutrino flux d(f) u (r') / dE u is given by Eq. (12) 
with Re being replaced by the distance from the Sun to 
the Earth, Rse- We use the Sun density profile given by 



p(r') = 236.93 g/cm 3 x exp(-10.098— ), (15) 

Rs 

where Rs is the radius of the Sun and we sum up all 
Sr' contributions until neutrinos reach the surface of the 
Sun. 



B. Upward Events and Muon Energy Loss 

High energy muons produced in neutrino charged- 
current interactions lose energy before they reach the de- 
tector as they travel through the rock or ice. The average 
energy loss of the muons with energy E over a distance 
dz during their passage through a medium with density 
p is given by, 



dE 

{^) = -{a + (3E)p, 



(16) 



where, a ~ 2 x 10~ 3 GeVcm 2 /g accounts for the ioniza- 
tion energy loss and (3 ~ 3.0 x 10~ 6 cm 2 /g accounts for 
the bremsstrahlung, pair production and photonuclear 
interactions. The parameter a is relatively insensitive to 
the composition of the material. The quantity [3 depends 
on composition of the medium and varies slowly with en- 
ergy [i^, [H, [4?| , and the average energy loss formula is 
not strictly applicable because of stochastic energy losses 
(47l [48| . For our purposes here, given the other uncer- 
tainties, using a constant (3 and approximating 



dE 

dz 



f dE 
"dz 1 



(17) 



With this assumption the initial energy at z = 0, E % , 
is related to the final energy E^ after traveling a distance 
z by 



E<Jz) = ef*"Ef + {e?»-l) 



(18) 



At low energies, for E n < 200 GeV, the contribution 
from P term is small (about 10 — 20%) and in this energy 
range, 

Ei(z) ~ El + apz . (19) 

Muons with energies of a few 100 GeV are stopped in the 
rock (p ~ 2.6 g/cm 3 ) before they decay. As an example, 
the stopping distance for 500 GeV muons is roughly 1 
km whereas the decay length, 7CT, for these muons turns 
out to be about 3000 km. For 50 GeV muons, the decay 
length is about 300 km, compared to a stopping distance 
of 100 m. 

The decay length information can still be included in 
the calculation by introducing the survival probability as 
the solution to the equation, 



dP m 



P., 



dE u E tl crp(a + PE^/rriy, 



(20) 



This leads us to the survival probability for a muon with 
initial energy E l and final energy E?, 



P SUIV (ELEf) 



'E[Y ( a + /3Ei\ r 



Ej, 



a + j3El 



(21) 



where T = m^/ (crap). 

With a distinction made between the energy of the 
muon when it is produced and the energy of the muon 
when it arrives at the detector for upward events, the 
formula for the upward muon flux is more complicated 
than Eq. (|T2"]) . Instead, we have 



dE u 



dr 



dE„ 



dP, 



cc 



dr dE], 



x dE~ v suiv[ ^ds; 

+ {v~>v) . (22) 

where E^ = E*. The minimum neutrino energy in the 
integral is E™ m = E^z), where z = R — r. The maxi- 
mum distance that muon travels to the detector and ends 
up with the final energy is R — R m in an d the relationship 
between E^ and E 1 ^ given by Eq. (16). 

Consider as a specific example annihilation at the core 
of the Earth. The detector is at a distance R = Re — 
6400 km from the core. In principle, Rmin = 0, how- 
ever, only muons produced near the detector will have 
sufficient energies to make it to the detector with an en- 
ergy above the detector threshold energy Eth- The muon 
average range is 



is sufficient. 
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Jp 



hi 



a + PEi 
a + (3E, 



th 



(23) 



5 



following from Eq. (| 18|) . For an initial muon energy of 
1 TeV, the muon average range is 1 km for a muon 
threshold energy of 50 GeV. In practice, then, R m in = 
Re — A, where A is, in general, less than 1 kilometer for 
the energies considered here. The upper limit for muon 
range is obtained by setting E 1 ^ — m x , in which case 
A = i? M (m x ,£: M ). 

A change of variable from r to z — Re ~ r yields a 
more familiar form of the integral for the muon flux from 
DM annihilations in the Earth's core, 



dE„ 



47Ti?| ,„ 



dzeP pz I dE, 



\dE v/ F ji 



x E» a + PE> \ (dot 



(24) 



Throughout Eq. (124|) , the initial muon energy is implic- 
itly a function of the final muon energy and the distance 
traveled, = E^{E^, z) from Eq. ([15]). 

Muon flux from DM annihilation in the Sun is given 

by 



dE„ 4ttR 2 se Jo 



dzePr* dE,/ 1 



dE, y 



E^ a + 0EI 
Eia + 0E„ 



TT exp(-p(r')<j C c$r' /m H )—^(r') 
L1 dE v 



(25) 



IV. RESULTS 

A. DM annihilation in the Earth's core 

To illustrate the muon flux's dependence on muon en- 
ergy, we begin with DM annihilation in the Earth's core. 
In addition to making a choice for m x , one must also 
make some assumptions about the cross section and main 
channel to produce neutrinos. For all of the figures for 
DM annihilation, we use a\ ~ 10~ 8 A^ pb [H| and the 
standard composition of the Earth as reviewed in Rcf. 

0- 

The upper curves in Fig. [T]show our results for XX ~ * 
v^Vp with B v = 1, for upward events (dot-dashed curve) 
and contained events (dashed curve) for m x — 500 GeV. 
The lower dot-dot-dashed and dot-dash-dashed curves in 
Fig. [T]come from XX t + t~ with B T = 1, followed by 
t — * VrfiUfi according to the energy distribution in the 
Appendix. We choose the tau channel as representative 
of the three body decays that also occurs in heavy flavor 
semileptonic decays. 



For direct production of neutrinos, d(j) v /dE v cx 5{m x — 
E v ). The cross section for neutrino production of muons 
smears the distribution. For contained events, one sees 
the smeared distribution directly in Fig. [1] Upward 
events have the additional energy redistribution from 
muon energy loss in transit that shifts the muon energy 
distribution to lower energies, enhancing the lower en- 
ergy flux relative to the contained flux, despite the fact 
that the range is shorter than 1=1 km. In the cascade of 
t — * Up — ► fi, shown with the lower curves, there is never 
a high energy peak and the upward events are always 



below the contained events for this value of 



Only 



for m x sufficiently higher than 1 TeV could the upward 
events be enhanced relative to the contained events. 

As an indication of the atmospheric neutrino back- 
ground, we also show upward (solid curve) and contained 
events (dotted curve) from a solid angle defined by a cone 
of half-angle 1° around the upward vertical direction. We 
use a simple parametrization for the flux of atmospheric 



■Vp (in units of GeV 1 km 

NoE^- 1 



-2-1-1 



sr ) from Ref . |4£ 



dE u dn 



1 + bE v cosQ 1 + eE v cost 



(26) 



where the values of the parameters iVo, 7 a, b, c and e are 
given in Table 1. The approximate angular resolution of 
the IceCube detector is 9 = 1°, however, DM annihilation 
can occur at angles larger than 1°. In particular, for a 
500 GeV neutralino, it has been shown [22| that most of 
the annihilation occurs within an angle of 9 ~ 2.7°. With 
this larger nadir angle, the solid angle for the atmospheric 
background is increased by a factor of ~ 8. 

The shape of the background atmospheric flux is very 
different from the signal of contained events for direct 
annihilation of DM into neutrinos. With the atmospheric 
contained events in the figure multiplied by a factor of 8, 
the contained event rate would dominate the background 
only for the high energy peak. Our sample calculation 
is for £>„(7q = 10 _8 A^ pb -1 , in which the capture and 
annihilation rates are in equilibrium. For the secondary 
neutrino production, from r decay, one needs B T <jf ~ 
10~ 7 A^ 4 pb -1 for this channel to be comparable to the 
background. Clearly measurements of the shape of the 
muon flux, both contained and upward, would be useful 
in searching for the dark matter signal. 

These values of the cross section a\ required for sig- 
nals on the order of the atmospheric background are 
sufficient for the condition for the equilibrium between 
capture and annihilation in the Earth's core to be sat- 
isfied. Even though only with significant enhancements 
of the capture rate (i.e. WIMP-nucleon cross section) 
or DM annihilation rate, the Earth might be a source of 
measurable rates for DM annihilation to neutrinos, it is 
a useful demonstration of the energy dependence of the 
muon flux. 

The energy dependence of Fig. [T] is at odds with the 
muon energy dependence sometimes found in the litera- 
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1.74 


a 


0.018 


b 


0.024 GeV" 1 


c 


0.0069 


e 


0.00139 GeV" 1 


No 


1.95 x 10 17 for v 
1.35 x 10 lr for V. 



TABLE I: Parameters for the atmospheric + f M flux, in 
units of GeV _1 km _2 yr _1 sr _1 . 
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-■ 500 GeV 

■ XX->w (upward) 
~ XX" >V V (contained) 
-■ XX->T + * ~ (upward) 

- ■ XX" >T+,C (contained) - 

■ ATM (contained) 

— ATM (upward) 




FIG. 1: Muon flux obtained from dark matter annihilation 
into neutrinos in the core of the Earth, when muons are 
created in neutrino interactions with nucleons in the rock 
below the detector (dot-dashed and dot-dot-dashed curves), 
when muons are created in the detector, i.e. contained events 
(dashed and dot-dash-dashed curves). The upper curves are 
for the direct production of neutrinos, while the lower curves 
are for neutrinos from tau decays. The background from con- 
tained atmospheric neutrinos, evaluated for a cone of angle 
9 = 1° are shown with the dotted (black) curve and the up- 
ward muon flux from atmospheric neutrinos is shown by the 
solid (black) curve. 



ture [id llll. l2ll. |41|. There, the upward flux of muons is 
written as 



dE u 



T A 
4ttR 2 



dE v 



E JE, 



dN v 



R^yE^, E t h) 



F 4 < 



(27) 



where Eth — 50 GeV. This expression accounts for the 
fact that muons have a range with an energy dependence, 
however, it does not account for the fact that over the 
distance i?(£' Al , E t h), the muon has a final energy of E t h- 
Eq. (|27p does not represent the energy dependent muon 
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pq 0.001 r 



m =500 GeV 

x i 

XX->w (Eq. (27)) : 
._. j^-ycV (Eq. (27)) ' 
— XX->tV(Ref. [22]) I 



300 400 500 

E (GeV) 



FIG. 2: Upward muons flux obtained using Eq. (27) for 
XX ~ * VpMfi (dashed curve) and for XX — * r+ r _ , followed 
by t — > Urfii/fi, (dot-dashed curve), and the muon upward flux 
for XX — * t + t~ channel from Ref. (solid curve). The up- 
ward muon flux from Eq. (27) is inconsistent with the upward 
flux shown in Fig. 1. 



flux, however, the integral number of upward events with 
Efi > E t h obtained using Eq. ([27]) and the results using 
Eq. (22) are approximately equal. In Fig. [2 we show 
the upward muon fluxes from Eq. (|27j) . for the direct 
neutrino production (dashed curve) and from the r de- 
cay (dot-dashed curve). Comparing results from Figs. 1 
and 2, we find that the upward muon flux of Eq. (1271) 
for XX ~> v v case follows more closely the contained 
muon flux at high energies presented in Fig.l (dashed 
curve) than the upward flux, with an enhancement at 
high because the muon range increases with muon 
energy. Clearly, the upward muon flux in Fig. 2 (dashed 
curve) does not accurately reflect the muon energy dis- 
tribution of upward events from DM annihilation in the 
Earth. Similarly, a comparison of upward muon flux for 
XX ~ * t + t", followed by r — > iv/ii^, obtained using 
Eq. (|27|) has a very different shape than the same flux 
obtained with Eq. (23). Comparable discrepancies are 
found between upward fluxes from Eq. (|27|) and our eval- 
uation of upward events for DM annihilation in the Sun 
as well. 

We also show with the solid line in Fig. [5] the re- 
sults for upward muon flux from the XX ~ * t + t~ from 
Ref. (22| . In Ref. [22( 1 . the flux of muons comes from 
a PYTHIA simulation of the resultant muon neutrino 
flux and a simulation of muon electromagnetic energy 
loss. A dark matter distribution has been assumed in the 
Earth's core and contribution from dark matter annihi- 
lation around the center of the core with specific angular 
cuts {9 < 5°) have been applied, so the normalization 
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should be lower. The energy distribution has qualita- 
tively the same behavior as our results, however, it does 
not vanish at the kinematic limit when = m x . 



B. DM annihilation in the Sun 

Similar conclusions can be derived in the case of cap- 
ture of WIMPs in the core of the Sun. As noted ear- 
lier, there is attenuation of the initial neutrino flux as 
it propagates from the core to the exterior of the Sun. 
The interaction length of the neutrinos with energy ~ 30 
GeV becomes equal to the column depth of the Sun (the 
average density of the core of the Sun is ~ 150 g/cm 3 ). 
At higher energies, the interaction length becomes even 
smaller and the neutrino flux is reduced significantly. We 
do not include neutrino oscillation in the Sun , which 
depending on the dark matter model, might affect the 
flux of + v^. 

In Fig. 3, we show the upward muon and the con- 
tained muon fluxes for the direct production and for the 
r production channels. In our calculations, we approxi- 
mate neutrino attenuation in the Sun with an exponen- 
tial suppression as presented in the previous section. We 
note that this effect becomes stronger for higher neutrino 
energies which manifests itself when m x is large. Recall 
that the charged current neutrino nucleon cross section 
increases with the neutrino energy. As an example, the 
muon flux decreases by a factor of 3 for m x = 250 GeV, 
factor of 10 for m x = 500 GeV and two orders of magni- 
tude for m x = 1 TeV, as compared to the case with no 
attenuation. 

We compare our results for muon flux with those in 
Ref. [12] , where there is assumption of dark matter dis- 
tribution in the core of the Sun and contribution from 
dark matter annihilation around the center of the core 
with specific angular cuts have been applied. Effects due 
to neutrino flavor oscillations in the Sun have not been 
incorporated. The shape of the energy distribution is 
similar to our result, but with lower normalization and 
with a lack of the kinematic cutoff when = m x . 

As in the case of the Earth, the upward muon flux from 
XX ~~ * y v is larger than the contained flux for muon en- 
ergies, Efj, < 380 GeV, while in the case when neutrinos 
are produced via XX ~ * t + t~, followed by r — > v T \iVn, 
the contained muon flux is always larger than the up- 
ward flux. We also show the angle-averaged atmospheric 
flux for a cone of half- angle 1°. For direct annihilation 
into neutrinos for the model in which the branching frac- 
tion is of the order of one, the signal is larger than the 
atmospheric background for both contained and upward 
muons. For the tau channel, signal is comparable to the 
background for upward muons when muons have energy 
around 200 GeV, however taking into account the effects 
of kinematics on the angular pointing of the muons at 
low energy may make this less apparent. 
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FIG. 3: Muon fluxes obtained from dark matter annihilation 
into neutrinos in the core of the Sun, for upward events (dot- 
dashed and dot-dot-dashed curves), and for contained events 
(dashed and dot-dash-dashed curves). The upper curves are 
for the direct production of neutrinos, while the lower curves 
are for neutrinos from tau decays. Background upward muons 
are shown with the solid (black) curve and the contained 
muons are shown with the dotted (black) curve, where the 
evaluation used the angle-averaged atmospheric neutrino flux 
integrated over a solid angle with 8 = 1°. The grey solid curve 
is from Edsjo's parameterization of the muon flux [221 ]. 



C. Muons in IceCube 

With the upward muon fluxes evaluated above from 
annihilation of DM in the Earth and the Sun, it is possi- 
ble to estimate the event rate of muons in IceCube using 
the muon effective area [43j]. Following Ref. [43|], we 
parameterized 

A eff (E^,9) ~ 2ttA (^)(0.92- 0.45 cos ff) (28) 

where 6 is the zenith angle measured from vertical and 

< 10 1 ' 6 GeV : 

Ao(E^) = 
10 1 ' 6 GeV < < 10 2 - 8 GeV : 

A (En) = 0.748(log 10 (^/GeV) - 1.6) km 

lO 2 ' 8 GeV < En : 

A (En) = 0.9 + 0.54(log 10 (£ M /GeV) -2.8) km . 

This effective muon area models the threshold detection 
effects near En ~ 50 GeV and local rock and ice below 
the IceCube detector [43T ]. 

To facilitate comparisons with other muon energy dis- 
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tributions which appear in the literature, we evaluate 



dE„ 



dE„ 



(29) 



for DM annihilation to neutrinos in the Earth and Sun 
which convert to muons outside the detector. Here, 
(A e ff) is the angle averaged effective area, averaged over 
zenith angles 9 — tt/2 — it. Fig. [4] shows our results 
for the upward muon flux times effective area with the 
solid and dot-dashed lines (solid for the Earth), and by 
comparison, the results for the contained muon flux mul- 
tiplied by 1 km 2 (dotted and dashed lines, dotted for 
the Earth) . The energy dependence of the effective area 
changes the shapes of the curves for upward muons at 
low energies, but it does not change the large discrepan- 
cies between our upward muon rates compared with Eq. 
(|2T|) at energies closer to E^ ~ m x . 



xx -> w 



> 
o 



m =500GeV 
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(Sun, upward) 
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FIG. 4: The upward muon flux times muon effective area 
obtained from dark matter annihilation to neutrinos in the 
core of the Earth (solid line) and Sun (dot-dashed line). For 
comparison, we also show the contained muon flux times 1 
km 2 for the Earth (dotted) and for the Sun (dashed). 



D. Cosmic Diffuse Neutrino Flux 



of halos of a given mass at a given redshift 42, 50, 51]. 
In Ref. [42| . dark matter annihilation process, XX ~ * vV, 
is proposed to be used to determine an upper limit on 
the annihilation cross section. 

The cosmic diffuse neutrinos for the XX ~ * v v channel 
from Ref. [421 ] is approximately a power law function of 
E„, i.e. 



/ d(j> v 
V dE v dn 



0.5 



A 



(£„/GeV) 
(m x /GeV) 3 - 5 



E v <m x . (30) 



In Ref. [42j], the normalization A is determined by set- 
ting the number of neutrinos from the diffuse flux (here 
approximated by Eq. (|30[) ) equal to the number of at- 
mospheric neutrinos from the same energy interval, from 
10~°- 5 m Y to m v , i-e. 



0.5 



dE v A 



v'lu 



(Ey/GcV) 
(m x /GeV) 3 - 5 



dE v 



dE„dtt 



where 



( d<t>„ \ 

\ dE v di1 I 



(31) 

dE i is the angle-averaged atmospheric flux 

v u / av 

given by Eq. (26). We consider here a case when m x = 1 
TeV. 

In Fig. we show upward and contained muon fluxes 
from cosmic diffuse neutrinos, integrating over the full 
2-7T solid angle. For E^ > 400 GeV, the contained flux 
dominates the upward flux, an artifact of the triangular 
shape of the neutrino flux. For comparison we also show 
the contained and upward fluxes for the atmospheric neu- 
trino background. The falling energy spectrum of the at- 
mospheric neutrinos results in the upward flux of muons 
from atmospheric neutrinos dominating the contained 
muon flux for > 400 GeV for m x = 1 TeV. 

The direct production of neutrinos XX ~* vP is t ne 
most favorable channel in terms of neutrino detection for 
the diffuse DM limits since the muon flux stands out more 
from the background than the muons from a x — > t — > ^ju 
cascade (or similar production and decay process) . In ad- 
dition, the XX ~ > y v channel has no other astrophysical 
observable. Nevertheless, the muon flux is not as dra- 
matic a peak in the falling neutrino induced atmospheric 
muon flux as the direct comparison of the neutrino fluxes 
is. A more comprehensive analysis of a diffuse DM an- 
nihilation signal could include both the vv and cascade 
channels as possibilities, and focus on the muon signals 
rather than the neutrino signals. 



In addition to the astrophysical object such as the Sun 
and the Earth being potential sources of dark matter, 
relic dark matter can also annihilate in halos in the uni- 
verse 42] . providing a promising source of cosmic diffuse 
neutrinos from dark matter annihilation. 

To determine this flux one needs to sum over all ha- 
los to yield a flux of neutrinos. This diffuse neutrino flux 
depends on several factors such as the evolution with red- 
shift, the radial density profiles and the number density 



V. CONCLUSIONS 

We have calculated muon fluxes from dark matter an- 
nihilation, when dark matter is trapped in the the Sun's 
(Earth's) core and when dark matter annihilates in ha- 
los in the universe (cosmic diffuse flux). Without us- 
ing a specific model for dark matter, we have considered 
XX ~ * y v an d XX ~ * t + t~, followed by r — > Vrfiv^ chan- 
nels as representatives of direct and of the secondary neu- 
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FIG. 5: Muon fluxes obtained from dark matter annihila- 
tion in the halos producing cosmic diffuse neutrinos: upward 
muons flux (dashed curve) and contained muon flux (dot- 
ted curve) compared with contained (dot-dashed curve) and 
upward (solid curve) muon fluxes from angle- averaged atmo- 
spheric neutrinos. Here, we take m x = 1 TeV. 



trino production. We have taken into account neutrino 
attenuation as it propagates from the core of the Sun to 
its surface. In the evaluation of the upward muon flux, 
we have incorporated muon energy loss, as described by 
the muon range. 

We have shown that our results exhibit a very different 
energy dependence than those obtained from Eq. ([77]) 
that is widely used in the literature PH HH, EU , however, 
there is reasonably good agreement with the parameter- 
ization of Ref. [22| away from the region of maximum 
energy ~ m x . Our results are obtained with the 
assumption that the dark matter annihilation occurs at 
the maximum rate, when the annihilation rate is half the 
capture rate. This is reasonable for the Sun but requires 
significant enhancement of the capture rate (or annihi- 
lation cross section) for the Earth to be in equilibrium 

In our calculation we used spin independent WIMP- 
nucleon cross sections which have much stronger exper- 
imental bound than the spin dependent cross sections 
[52| . In the core of the Sun the capture rate might be 
dominated by the spin dependent (SD) WIMP-hydrogen 
nuclei interactions, which would increase the signal rates 
by a couple of orders of magnitude and still be consis- 
tent with Amanda limits on annihilation rates [53| . In 
the dark matter model in which there is a low velocity 
enhancement of the DM annihilation cross section [271 ] . 
introduced as an explanation for the positron excess ob- 
served in cosmic ray experiments [30, HH, HH, |33{ . it is 
possible for the WIMPs in the core of the Earth to be in 



the equilibrium as well. 

Furthermore, incorporating neutrino oscillations and 
the regeneration effects in the Sun will likely affect the 
final muon flux especially in the models which possess 
an asymmetry in the initial neutrino fluxes or where 
XX ~ * t + t~ is the dominant mode [l(J. We have used 
a model independent normalization, <7q ~ 10 _8 A^ pb 
and Bp = 1 to evaluate the muon flux. We find that 
for this branching fraction signals from XX ~> v ^ an d 
XX ~ * T+T ~j followed by r — > Vt^lv^, when DM annihi- 
lation happens in the core of the Sun, are comparable or 
even larger than the background (upward) muons from 
atmospheric neutrinos. In the case of direct neutrino 
production, the upward muon flux is larger than the con- 
tained flux for Ef, < 350 GeV for m x = 500 GeV, due 
to the muon range. When neutrinos are produced via 
secondary processes, contained events always dominate 
upward muons. 

Cosmic diffuse neutrinos, produced directly in DM an- 
nihilation, give a very weak energy dependence of the 
contained and upward muon flux, in contrast to the steep 
energy dependence of the atmospheric background. The 
upward muon flux from cosmic diffuse neutrinos is dom- 
inant over the contained flux for muon energies below 
400 GeV for m x — 1 TeV, which also happens to be the 
energy at which the signal becomes dominant over the 
background muons from atmospheric neutrinos. 

Model dependence is an important element, for exam- 
ple, XX ~ * v v is n °t allowed for DM at rest when the 
DM particles are neutralinos [l(| ■ However, with the 
formalism developed, one can determine muon fluxes for 
specific dark matter model by summing up the contribu- 
tions from all decay channels weighted with correspond- 
ing branching fractions 54]. Thus, measurements of the 
muon energy distribution in neutrino telescopes, such as 
IceCUBE and KM3, could provide valuable information 
about the origin of the dark matter sector and funda- 
mental properties such as the dark matter mass and its 
couplings. 
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VI. APPENDIX:MUON NEUTRINO 
DISTRIBUTIONS 

A. Neutrino energy distribution from direct 
production: 



to be the same. The decay branching fraction is denoted 
by Bf for a given decay mode f,f = v,r,b,c. The b and 
c quarks hadronize before they decay into neutrinos. The 
hadronization effect is taken into account by scaling the 
initial quark energy, Ei n , in the form Ed — ZfEm, where 
Zf = 0.73, 0.58 for b and c quarks, respectively [551 ]. 

Neutrino energy distribution from the decay of 
t + t~, bb, cc is approximately 



Neutrino energy distribution when neutrinos are pro- 
duced directly from dark matter annihilation is given by, 



dN„ 
dE v 



S(E V — m x ) 



(32) 



B. Neutrino energy distribution from t + t ,bb,cc 
decay modes: 



dN v 
dE v 

where 



2B 



{E in , B f ) 



2x 6 



where x 



Ej T 



< 1 



(33) 



(m x , 0.18) r decay, 
(0.73m x , 0.103) b decay, (34) 
(0.58to x , 0.13) c decay. 



In these decay modes, we use the unpolarized decay 
distributions, so the v and V distributions are assumed 
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